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We investigate the effective interaction between two randomly charged but otherwise net-neutral, 
planar dielectric slabs immersed in an asymmetric Coulomb fluid containing a mixture of mobile 
monovalent and multivalent ions. The presence of charge disorder on the apposed bounding sur¬ 
faces of the slabs leads to substantial qualitative changes in the way they interact, as compared with 
the standard picture provided by the van der Waals and image-induced, ion-depletion interactions. 

While, the latter predict purely attractive interactions between strictly neutral slabs, we show that 
the combined effects from surface charge disorder, image depletion, Debye (or salt) screening and 
also, in particular, their coupling with multivalent ions, give rise to a more diverse behavior for 
the effective interaction between net-neutral slabs. Disorder effects show large variation depending 
on the properly quantified strength of disorder, leading either to non-monotonic effective interac¬ 
tion with both repulsive and attractive branches when the surface charges are weakly disordered 
(small disorder variance) or to a dominating attractive interaction that is larger both in its range 
and magnitude than what is predicted from the van der Waals and image-induced, ion-depletion 
interactions, when the surfaces are strongly disordered (large disorder variance). 


I. INTRODUCTION 

Interactions between neutral dielectric bodies are tra¬ 
ditionally viewed as being due predominantly to electro¬ 
magnetic field fluctuations, or equivalently, dipole fluctu¬ 
ations that give rise to van der Waals (vdW) fluctuation- 
induced interactions between them [l|, [^. These inter¬ 
actions are attractive between identical bodies in vac¬ 
uum or in a polarizable medium, such as water or 
an aqueous electrolyte (or Coulomb fluid). They con¬ 
tribute one of the two key ingredients of the classical 
Derjaguin-Landau-Verwey-Overbeek (DLVO) theory of 
colloidal stability, the other one being the mean-field elec¬ 
trostatic interactioUjwhich is repulsive for like-charged 
colloidal surfaces [J-Q- 

Recent works have, however, highlighted the role of 
image-induced, ion-depletion effects in this scenario (see, 
e.g.. Refs, [sl-l^ and references therein), leading to de¬ 
pletion of mobile solution ions from the vicinity of dielec¬ 
tric interfaces and, therefore, to an additional attractive 
force between apposed dielectric boundaries. This is be¬ 
cause most dielectric surfaces in the context of bio- and 
soft materials have a lower (static) dielectric constant 
than that of water and, therefore, solution ions expe¬ 
rience repulsion from their same-sign image charges in 
the proximity of dielectric boundaries [s^. The recent 
advances in the study of image-induced, ion-depletion ef¬ 
fects follow on the trail of the Onsager-Samaras frame¬ 
work formulated originally in the context of the surface 
tension problem of electrolytes [i^, |4l| . Such non-mean- 
field effects, which belong to the general class of depletion 
interactions [d^, arise due to the discrete nature of mo¬ 
bile ions neglected in the collective mean-field description 


* Email: a.najiQipm.ir 


based on the standard Poisson-Boltzmann theory 
The studies of image-induced, ion-depletion effects have 
been focused exclusively on the case of strictly neutral 
(charge-free) dielectric surfaces. In this case, the ion- 
depletion interactions can be amplified in the presence 
of mobile multivalent ions in the solution [l^, [l^ due to 
stronger ion-image repulsions for these ions, even when 
the multivalent ions are present at small bulk concentra¬ 
tions around just a few mM 

In this paper, we revisit the problem of interaction be¬ 
tween neutral dielectric surfaces in a Coulomb fluid by 
adding to it a novel feature: We relax the constraint of 
strict electroneutrality of surface boundaries, considered 
so far in the literature, by assuming that the surfaces are 
neutral only on the average, while microscopically they 
carry a quenched (fixed) random distribution of positive 
and negative charges. We show that this seemingly sim¬ 
ple generalization leads to significant qualitative changes 
in the distribution of ions and, consequently, also in the 
effective interactions between dielectric surfaces, espe¬ 
cially when the intervening Coulomb fluid contains mo¬ 
bile multivalent ions. 

Disordered charged systems are abundant in soft mat¬ 
ter with examples ranging from polycry stalline surfaces 
with patchy surface potentials [4^l445l| . dielectric con¬ 
tact surfaces l46|. y apor-deposited amorphous films on 
solid substrates l47ll. s urfactant-coated surfaces [isl - l^ . 
DNA microarravs |53l. [s^, intrinsically disordered pro¬ 
teins [H, [5^, patchy colloids and random poly¬ 
electrolytes and polyampholytes |58l - l60l| . Surface charge 
disorder in these examples can exhibit highly random 
distributions as well as patchy and heterogeneous pat¬ 
terns, originating from different sources including spe¬ 
cific electronic and/or structural properties of materi¬ 
als involved, surface grafting or adsorption of charged 
molecules and/or contaminants, synthetic and fabrica- 



2 


tion processes, etc. The surface charge disorder can be 
highly sample specific and, at the same time, can de¬ 
pend strongly on the method of preparation. It may be 
set and quenched for each sample (which is the case of 
interest in this paper), annealed (in which case the sur¬ 
face charges are mobile and in thermal equilibrium with 
the rest of the system), or partially quenched or partially 
annealed [^,[62. Motivated by these examples, study 
of charge disorder and, in particular, effective interac¬ 
tions between random charge distributions has witnessed 
growing attention from theoreticians over the last sev¬ 
eral years [H, Issl - fs^l , as well as from the simulation side 
where initial steps have been taken to include the effects 
of charge heterogeneity and disorder [HtI - Ib^ . It^ [85l - l^ . 
Disorder effects have been studied extensively in the con¬ 
text of electromagnetic fluctuation-induced interactions 
between two app osed, randomly charged surfaces in vac¬ 
uum [3, [7^82l |. where they have been associated with 
anomalously long-ranged surface interactions observed in 
recent experiments [dj, I3 • They have also been studied 
in situations where a weakly coupled Coulomb fluid, con¬ 
taining, e.g., monovalent cation and anions, intervenes 
between the bounding surfaces [H, HH, IbSl - f^ It^ . 
This however leaves out the case of asymmetric Coulomb 
fluids containing both monovalent ions and multivalent 
(counter-) ions. These kinds of systems are in fact quite 
common in experiments within the biological context as 
in the case of viruses, DNA condensates or other charged 
biopolymer aggregates [ssl - l^ and are expected to be¬ 
have very differently since multivalent ions are known to 
couple strongly with fixed surface charges. 


Strong-coupling behavior of multivalent counterions 
at uniformly charged surfaces or surfaces with regular 
charge patterns has been studied extensively over the 
last decade and its connection to exotic phenomena such 
as like-charged attraction has been throughly discussed 
(for recent reviews and a more exhaustive list of refer¬ 
ences, see Refs. BmiiMoi). In disordered charged 
systems, such strong-coupling phenomena have been con¬ 
sidered only in a few cases so far [tI, HI, by as¬ 
suming that bounding surfaces carry a finite mean surface 
charge density, to which multivalent counterions can cou¬ 
ple strongly, in the same sense as considered in the con¬ 
text of uniformly charged surfaces as noted above. Yet, 
the presence of charge randomness on bounding surfaces 
was shown to give rise to novel phenomena such as strong 
surface attraction of multivalent counterions, character¬ 
ized by a density profile that diverges at the surface, in 
clear violation of the contact-value theorem established 
for uniformly charged surfaces. This kind of behavior 
originates from a singular, attractive, single-ion poten¬ 
tial, which is created by the randomness in the distribu¬ 
tion of surface charges and, as such, depends on the sur¬ 
face charge variance. Consequently, one can also show 
that the thermal entropy of counterions is diminished 
upon introducing a finite degree of charge randomness 
on the boundaries. In other words, the system becomes 
more ‘ordered’ as a direct outcome of the interplay be¬ 


tween the thermal entropy of ions and the configurational 
entropy of charge randomness, entering through an en¬ 
semble average over various realizations of the disordered 
charge distribution. This peculiar disorder-induced ef¬ 
fect, which stands at odds with what one m ay ex pect 
intuitively, has been referred to as antifragility [l05j| . 

Our analysis in this paper is focused on yet another 
facet of the disorder-induced effects by assuming, in a 
system of two plane-parallel dielectric slabs, that the ran¬ 
domly charged surfaces of the slabs bear no net charge. 
This eliminates the direct strong coupling between mul¬ 
tivalent ions and the mean surface charge and, thereby, 
also the ensuing stron g-coupl ing interactions, considered 
in our previous papers [83, IS^, that would otherwise com¬ 
pletely mask the vdW and image-induced, ion-depletion 
interactions between the slabs in an asymmetric Coulomb 
fluid. This allows us to address the question of how the 
presence of surface charge disorder affects the standard 
picture for the interaction of neutral bodies based on the 
vdW and ion-depletion effects. 

The surface charge disorder has several different impli¬ 
cations: First, it contributes a repulsive interaction be¬ 
tween the slabs, which comes from self-interactions of dis¬ 
order charges on the bounding surfaces with their image 
charges; this contribution counteracts the attractive vdW 
interaction as discussed thoroughly elsewhere [tsI - I^ . 
Then, as noted above, the individual ions in the Coulomb 
fluid experience an attractive disorder-induced potential, 
which strongly attracts them toward the bounding sur¬ 
faces. This effect counteracts the ion-image repulsions 
that tend to deplete ions from the slit region between 
the slabs and form the basis of the image-induced, ion- 
depletion mechanism for attraction between strictly neu¬ 
tral slabs. Thus, randomly charged bounding surfaces 
tend to accumulate more mono- and multivalent ions in 
the slit. The system, however, responds differently to 
the increased number of ions: While the osmotic pressure 
due to monovalent ions increases and even becomes re- 
pulsive, consistent with the standard mean-field picture, 
the osmotic pressure due to multivalent ions, becomes 
ever more attractivel This behavior is rooted in the com¬ 
bined effect of the surface charge disorder and the pres¬ 
ence of mobile multivalent ions, with the latter creating 
strong inter-surface attractions upon further accumula¬ 
tion in the slit. As a result, the effective total inter¬ 
action pressure between randomly charged, net-neutral 
dielectric surfaces can differ qualitatively from what one 
expects based on the standard vdW [ij, Q and image- 
induced, ion-depletion interactions in the case of strictly 
neutral surfaces (see, e.g., Refs. [§-[3^ and references 
therein). 

The net effect due to the interplay between disorder, 
mono- and multivalent-ion contributions yields a qualita¬ 
tively different behavior for the effective surface-surface 
interactions, depending on the strength of disorder as 
quantified b y th e disorder coupling (or strength) param¬ 
eter [^, ISMISJI ■ This feature bears some conceptual 
resemblance to the strong-weak coupling dichotomy that 
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FIG. 1. (Color online) Two infinite, plane-parallel dielectric 
slabs of dielectric constant tp and randomly charged inner 
surfaces are immersed in a bathing ionic solution of dielectric 
constant em- The solution contains a mixture of monovalent 
and multivalent salts. Multivalent ions are confined in the 
slit region and are shown by large spheres; monovalent salt 
anions and cations are shown by small red and blue spheres. 
The slabs are assumed to be neutral on the average and their 
thickness is taken to be infinite. 

exists for net-charged surfaces, dependen t on the elec¬ 
trostatic coupling parameter in that case |l0fil4l0^ . In 
weakly disordered systems, corresponding to a small dis¬ 
order coupling parameter, one then discerns a distinct 
non-monotonic behavior for the interaction pressure as 
a function of the separation between the slabs, with a 
pronounced repulsive hump at intermediate separations; 
conversely, in strongly disordered systems, correspond¬ 
ing to a large disorder coupling parameter, the effective 
interaction pressure becomes strongly attractive, with a 
range and magnitude larger than that of the vdW or 
the image-induced, ion-depletion interaction pressure as 
found between strictly neutral surfaces. 

The organization of the paper is as follows: In Section 
im we introduce our model and, in Section Hill we briefly 
discuss the theoretical background and present the gen¬ 
eral results for the two-slab system. The results of our 
analysis for the density profile of multivalent ions and the 
effective interactions between the slabs are presented in 
Section |TVl We conclude our discussion in Section El 


II. THE MODEL 

Our model consists of two plane-parallel dielectric slabs 
of infinite surface area S and dielectric constant Cp with 
inner surfaces placed perpendicular to the z axis at a 
separation distance of d (see Fig. [T]). The slabs are im¬ 
mersed in an asymmetric Coulomb fluid of dielectric con¬ 
stant Cm at room temperature T. The Coulomb fluid is 
a mixture of a monovalent I : I salt of bulk (reservoir) 
concentration rig and a multivalent q : 1 salt of bulk con¬ 
centration Co with multivalent ions having charge valency 
of g > 0. The multivalent ions are confined within the 
slit region —d/2 < z < d/2, while the monovalent ions 
are dispersed throughout the space except in the region 
occupied by the dielectric slabs that are assumed to be 


impermeable to mobile ions |l09l| . We are interested only 
in the cases where the slab thickness is much larger than 
the Debye (or salt) screening length and, thus, in the cal¬ 
culations to be presented later, we shall assume that the 
slab thickness is infinite. 

The inner surfaces of the slabs are assumed to bear 
a quenched, random charge distribution p(r), while they 
remain electrically net-neutral. We assume that the ran¬ 
dom surface charges are distributed according to a Gaus¬ 
sian probability distribution function, which is deter¬ 
mined fully by its two moments, ((p(r))) = 0, and 

((p(r)p(r'))) = g(r)S(r - r'). (1) 

with 

ff(r) = geo[d(z-l-d/2 )-l-S(z-d/2)], (2) 

where g > 0 is the surface disorder variance. By assump¬ 
tion, thus, the disordered charge distributions on the two 
slabs are statistically uncorrelated. The effects of sur¬ 
face charge correlations (or “patchiness”) and the inter¬ 
nal structure of multivalent ions, which are a ssumed t o 
be point-like here, will be considered elsewhere [na[mi. 

III. INTERACTION EREE ENERGY: GENERAL 
RESULTS 

A. Formal background 

In the present model, we have two types of mobile ions, 
namely, mono- and multivalent ions, that are expected 
to behave very differently in the presence of dielectric 
boundaries. On general grounds, and if the surfaces are 
assumed to be strictly neutral (charge free), mobile ions 
interact only with their image charges. This interaction 
scales with the second power of the ionic valency, be¬ 
ing thus much larger in the case of multivalent ions with 
q ^ 1 (e.g., trivalent and tetravalent ions). In this latter 
case, the ion-image interactions are dominated by the 
self-image interaction. Note that in most realistic ex¬ 
amples, ions are dissolved in water, which is a highly 
polarizable medium with = 80 at room temperature 
(T = 293 K) and, thus, we often have €m > Cp, giving 
image charges of the same sign and, hence, repulsive self- 
image interactions for each ion in the slit. Such repulsive 
interactions lead to statistical correlations (manifested as 
ion depletion) between the mobile ions and the bound¬ 
ing surfaces that are expected to be weak for monovalent 
ions, but quite strong for multivalent ions. This leads 
to a complicated problem in which different ionic species 
show distinctly different couplings to local electrostatic 
fields generated by the boundaries. While the mono¬ 
valent ions can be therefore described appropriately by 
mean-field-type theories, such as the Poisson-Boltzmann 
or the Debye-Hiickel theory , the multivalent ions 
require an altogether different description that should ac¬ 
count for such large correlation effects. 
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The situation is, therefore, analogous to the one found 
in the case of asymmetric Coulomb fluids confined be¬ 
tween charged surfaces, where the dominant factor deter¬ 
mining the behavior of multivalent counterions is again 
their electrostatic coupling to the local fields that are 
generated by the boundaries. The difference is however 
that in the the case of a charged surface the coupling 
is deter mined by the so-called electrostatic coupling pa¬ 
rameter flObUIo^ . depending on the charge density of 
the surface boundaries, enabling the so-called strong¬ 
coupling approximation for multivalent cou nterions a t 
charred surfa ces (see, e.g., Refs. 0, m IggHTol Hoi- 
[ioi,[iiMi3 and references therein ), o r mor e gen e rally , 
the dressed multivalent-ion theory [23, Eq! [ilO, [ml, 
which provides a very good approximation for the study 
of asymmetric Coulomb fluids over a wide range of 
parameters as verified by explicit-ion and implicit-ion 
simulations. The dressed multivalent-ion theory repro¬ 
duces the mean-field Debye-Hiickel theory and the stan- 
dar d strong-coupl i ng theory for counterion-only systems 
' 9 ^ 10.11 llOfil - llok Ill2l - lll8l | as two limiting theories in 


the regime of large and small Debye screening lengths 
and , therefore, br idges the gap between these two limits 
[13, il03l ll20Ml23j| . The key step in this latter approach 
is to integrate out the degrees of freedom associated with 
monovalent ions by means of a linearization approxima¬ 
tion, justified only for highly asymmetric Coulomb flu¬ 
ids q » 1, and yielding an effective Debye-Hiickel (DH) 
interaction between the rema ining multivalent ions and 
the surface charges (if any) |l2fll |. Then, one expands 
the partition function of the system in terms of the fu- 
gacity (or bulk concentration, Ac = cq) of multivalent 
ions or in terms of the inverse electro stati c coupling pa¬ 
rameter for counterion-only systems |lOfil |. The leading 
order of the virial expansion can be cast into a sim¬ 
ple analytic theory because of its single-particle struc¬ 
ture, which can successfully describe various features 
of strongly charged systems containing multivalent ions 
(see, e.g. Refs. and 

references therein). The regime of applicability of this 
theory has been discussed exten sively in recent literature 
0.[^[^[^ [Tfn[ - [IM[iM[^^ , which we therefore 

do not reiterate here. 

In the case of strictly neutral (charge-free) surfaces, 
this strategy was shown to be effective as well 27| , since, 
as is often the case in experimental systems containing 
asymmetric ionic mixtures [88l - [^ . multivalent counteri¬ 
ons of high valency (for instance, CoHex^+ or polyamines 
such as Spd^+, Sp'^’*') are present only in small bulk con¬ 
centrations, of about a few mM, justifying fully the virial 
expansion scheme that underlies the dressed multivalent- 
ion theory. This is the approach of choice that we adopt 
also in the present context, where the surface boundaries 
carry random charges and are neutral only on the aver¬ 
age. The only additional step here is that, due to the 
quenched disorder in the surface charge distribution, the 
free energy of the system follows by averaging over the 
whole ensemble of charge distributions, p(r) [sj, [13|- 


B. Planar slabs 


Without delving further in to details avai lable in pre¬ 
vious publications [13, [US EM EM Em, we proceed 
by giving the general expression for the free energy of the 
two-slab system considered in this work, which has the 
form [m 

J-= iTr[g(r)G(r,r')] - AefeT J dr D(r) (3) 

Here, /? = l/{kBT) and the Green’s function of the sys¬ 
tem is defined via 


-eoV-e(r)VG(r, r')-|-eoe(r)At^(r)G(r, r') = (5(r-r'), (4) 


where K(r) is the Debye (or salt) screening parameter, 
which is non-zero only outside the region occupied by 
the dielectric slabs, i.e., = dTr/ent, with Ub = 2nQ-\-qco 

being the total bulk concentration of monovalent ions. 

The first term in Eq. ([3]) represents the free energy of 
the system in the absence of multivalent ions. It stems 
from self-interactions of disorder charges on the bounding 
surfaces of the slabs with their image charges (note that 
the random charge distributions on the two slab are un¬ 
correlated and, on the average, do not interact with each 
other). This term depends only on the disorder variance, 
g, and is non-vanishing only in inhomogeneous systems 
with a finite dielectric discontinuity at the bounding sur¬ 
faces and/or a spatially inhomogeneous distribution of 
salt ions. This contribution has been analyzed in the 
context of the fluctuation-induced forces between disor¬ 
dered surfaces in vacuum or in a weakly coupled Coulomb 

fluid [zilill- 


The second term in Eq. ([S]) represents the contribu¬ 
tion from multivalent ions on the leading (virial) order, 
in which ufr) is the effective single-particle interaction 
energy [M, Ml 


2 2 2 2 
M(r) = ^^Gim(r,r) -/3^— 


dr'g(r')[G(r,r')]^. (5) 


The first term in the above expression represents the 
self-image interaction of individual multivalent ions in 
the slit and the second term represents the contribu¬ 
tion of the surface charge disorder to the effective single¬ 
particle interaction energy. In the first term, Giin(r,r) = 
G(r, r) —Go(r, r) gives the dielectric and/or salt polariza¬ 
tion, or the image-charges effects (corresponding to the 
generalized Born energy), in which the free-space Green’s 
function, Go (r, r) (corresponding to the formation en¬ 
ergy of individual ions in a homogeneous background), 
defined via —eoem(V^ — K^)Go(r,r') = (5(r — r'), is sub¬ 
tracted from the total Green’s function. The second term 
in Eq. (© is found to be proportional to the disorder vari¬ 
ance and shows an explicit temperature dependence and 
a quadratic dependence on the Green’s function and the 
multivalent-ion charge valency, q (these latter features 
can be understood by noting that the disorder term in¬ 
deed stems from the sample-to-sample fluctuations, or 
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variance, of the sample-dependent single-particle inter¬ 
action energy as discussed in Ref. [13 )■ 

In the strong-co upling lirnit or within the multivalent 
dressed-ion theory [l06l llOSL Il20l| , the number density of 
multivalent counterions can be obtained in terms of the 
effective single-particle interaction energy as [H, H^l 

c(r) = A,fI(r)e-^“«. (6) 


In the specific example of two planar slabs, we can take 
advantage of the translational invariance of the Green’s 
function with respect to transverse (in-plane) direction 
coordinates p = {x,y) and p' = {x\y'), since G(r,r') = 
G{p, p'; z, z') is only a function of \p — p'\, z and z\ 
and write the free energy in terms of its Fourier-Bessel 
transform G{Q;z,z') defined through 


G(r,r') = 


/o 


QdQ 

2tt 


G{Q-z,z')MQ\p-p'\). (7) 


For two semi-infinite slabs, one has 

W 

X (e^‘^cosh7(z -|- z') -t- Ag cosh7(2; — z'))]. 
where 7 ^ = -I- and 


Ag = 


e-ml + f^pQ 


(9) 


In the absence of salt screening {k = 0), Ag reduces to 
the bare dielectric discontinuity parameter 


A = (10) 

“r Cp 

which gives a measure of the magnitude of “dielectric- 
image” charges; while in a dielectrically homogeneous sys¬ 
tem, we have 


Ag 



1-Q 
7 + Q’ 


( 11 ) 


The interaction free energy of the system (per and 
per unit area, S) can be written as 


ojr j.d/2 

= glBfin,d,A) - dzc{z), 

* d-d/2 


(14) 


where we have subtracted additive terms that are inde¬ 
pendent of surface separation d and defined 


f{K,d,A) 



Ag(l + Ag)^ 

7(e2d7 _ A^)' 


(15) 


IV. RESULTS 

A. Distribution of multivalent ions 

Let us first focus on the case of a dielectrically homo¬ 
geneous system with A = 0. The density profiles of mul¬ 
tivalent ions for this case are shown in Fig. [5^, where we 
have rescaled the density profiles with their bulk value 
Co and the position in the slit —d/2 < z < d/2 with 
the Bjerrum length /b- The inter-surface distance and 
the screening parameter are rescaled in the same way. 
In the figure, we have fixed d/l^ = 5 and kZb = 0.35, 
equivalent to mono- and multivalent salt concentrations 
of no = 20 mM and cq = 1 mM when the Bjerrum length 
is taken as Ib = 0.71 nm (appropriate for water at room 
temperature, i.e., with T = 293 K and e™ = 80) and 
the multivalent ion valency is taken as q = 4 (note that 
throughout this paper we focus on the case of asymmetric 
Coulomb fluids with tetravalent ions but the generaliza¬ 
tion of our results to other values of q is straightforward). 

The surface charge disorder variance is shown in the 
figure in terms of the dimensionless disorder coupling (or 
strength) parameter [s^, [s^l 

X = 27rg^/|g, (16) 


which gives a measure of the salt-induced image effects, 
or loosely speaking, “salt-image” charges. Both image 
charge effects lead to depletion of ions from the vicinity 
of dielectric surfaces in a medium of relatively large polar¬ 
izability, but they exhibit some fundamental differences, 
which have been discussed in the context of charged sur¬ 
faces in our previous works [H, HJl . 

In the two-slab geometry, we can then re-express the 
(number) density profile of multivalent ions in the slit 
region, —d/2 < z < d/2, as 

c{z) = COB- (12) 

where 

22 rco 

u{z) = —-^ QdQG(Q;z,z) 

47r Jo 

J°°QdQ [G^ (Q; Z, -d/2) + G^Q- z, d/2)]. 

(13) 


which is varied in the figure as % = 0,1 and 2 (corre¬ 
sponding to strictly neutral surfaces with g = 0, and 
randomly charged surfaces with g = 0.02 nm“^ and 
0.04 nm“^, respectively). 

As seen, for strictly neutral (or disorder-free) di¬ 
electrics, the density profile (black solid curve) shows a 
partial depletion of multivalent ions from the vicinity of 
the surface boundaries. This is caused by the salt-image 
repulsion (corresponding to the first term in Eq. ([5|)), 
which, as noted in the previous Section, are caused by 
the discontinuity in the distribution of salt ions across 
the dielectric interfaces, giving a non-vanishing Ag ac¬ 
cording to Eq. (HB). However, when the surfaces are 
randomly charged (dashed curves), the situation turns 
out significantly different and the multivalent ions are 
attracted quite strongly towards the surfaces despite the 
“soft” salt-image repulsions (moreover, one can note from 
the shown profiles that a larger amount of multivalent 
ions are pulled into the slit from the bulk solution). In 
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FIG. 2. Rescaled density profile of dressed multivalent ions as a function of the rescaled normal position in the slit between the 
randomly charged inner surfaces of two net-neutral dielectric slabs for (a) dielectrically homogeneous system with A = 0 and 
(b) dielectrically inhomogeneous system with A = 0.95 (appropriate for the water/hydrocarbon interface). Other parameters 
are fixed as d/lB = 5 and kIb = 0.35 (cq = 1 mM, no = 20 mM). Panel (c) shows the threshold value, x* (above which the 
unimodal density profile changes to a bimodal one) as a function of A for d/ls = 5, co = 1 mM and no = 20, 40 and 60 mM. 


fact, the resulting disorder potential acting on individ¬ 
ual multivalent ions (corresponding to the second term in 
Eq. ([5])) exhibits a singular (logarithmic) behavior at the 
two surfaces and gives an algebraically diverging density 
as c(z) ^ (d^/4 — when ^ —>• rkdj^. This kind 

of phenomena has been discussed in detail in the con¬ 
text of disordered charged surfaces bearing a net charge 
density [13, , where we show that the singular behav¬ 

ior of the density profile is closely connected with the 
anti-fragile behavior of multivalent counterions; Adding 
a quenched disorder component to an otherwise uniform 
distribution of surface charges leads to a lowered entropy 
(or thermal disorder) for the counterions. In the present 
context, with net-neutral surfaces, the disorder-induced 
effects are qualitatively similar as they are regulated only 
by the variance of the charge distribution, g, rather than 
the net charge density of the surface. 

Finite dielectric discontinuity at the bounding surfaces, 
A > 0, creates a much stronger depletion effect than in 
the case of dielectrically homogeneous system (see Fig. 
[^b)- In fact, the density profile of multivalent ions van¬ 
ishes in the immediate vicinity of the surfaces because 
of the strong dielectric-image repulsion, provided again 
by the first term in Eq. (0. This term is singular itself 
and diverges on approach to the surface as ~ l/(d/2 =F z) 
when z —^ '^dl2, thus dominating over the singular po¬ 
tential created by the charge randomness on the bound¬ 
aries and, leading to the vanishing contact density of 
multivalent ions regardless of the disorder strength, y, 
as seen in the figure (even though a larger amount of 
multivalent ions are found in the slit at larger values of 
x). The competition between these two mechanisms of 
repulsion and attraction leads to a change in the shape 
of the density profile from unimodal to bimodal beyond 
a threshold value of the disorder strength, y*. The data 
in Fig. ^ show x* as a function of the dielectric discon¬ 
tinuity parameter. A, for fixed d/l^ = 5, cq = 1 mM and 


different values of uq = 20,40 and 60 mM (corresponding 
to kIq = 0.35,0.48 and 0.58, respectively). Clearly, for 
larger salt screening and/or dielectric discontinuity pa¬ 
rameter, a larger degree of charge disorder is required in 
order to counteract the image-charge repulsions. These 
features of the density profiles are qualitatively similar to 
those found in the case of charged surfaces [8^ and, 
therefore, we shall not delve further into the details of 
the behavior of the density profile and proceed with the 
analysis of the effective interaction between net-neutral 
surfaces. 


B. Effective interactions 


Effective interactions between neutral dielectric slabs 
are standardly described in terms of the vdW interactions 
as, for instance, formulated within the Lifshitz theory P, 
P|. The vdW interaction messure for two plane-parallel 
slabs can be expressed as [l[ 


Pvdw = —k'aT 


' QdQ yAg A 

27r 1 — A? dird^ 


(17) 


where the first term comes from the thermal zero- 
frequency mode of the electromagnetic field-fluctuations 
and the second term comes from the higher-order Mat- 
subara frequencies, pertaining to quantum fluctuations. 
A is then the quantum part of the Hamaker coefficient, 
which, based on an upper bound estimate in the case of 
hydrocarbon slabs interac ting acro ss an aqueous medium, 
can be taken as A = 3 zJ [27Lll24| . The vdW interaction 
is dominant mostly at the scal e of a few nanometers for 
the inter-surface distance [l04j . 

In the present model, one needs to account for the elec¬ 
trostatic contribution to the inter-surface pressure, Pes, 
as well. This contribution can be written as the difference 
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between the slit pressure of multivalent and monovalent 
ions and the bulk pressure Pb = {ub + co)kBT, where 
Ub = 2no + qcQ is the total bulk concentration of mono¬ 
valent ions, i.e., Peg = Pg — Pb- The slit pressure can be 
calculated as Pg = —dP/{Sdd) + n{0)kBT, in which the 
first term is the contribution of multivalent ions that fol¬ 
lows from the derivative of the free energy expression, Eq. 
(HI, w.r.t. the inter-surface distance when all other pa¬ 
rameters are kept constant, while the second term is the 
contribution of monovalent ions that has been expressed 
in terms of the total mid-plane density of monovalent ions 
27| . This latter quantity can be estimated here through 
the relation n{z) = rif, exp[—/3it(2:)] for z = 0, which 
has been shown by means of explicit Monte -Carlo sim¬ 
ulations of a disorder-free system [23, Il03j to give an 
accurate estimate of the distribution of monovalent ions 
in the slit. 

The electrostatic pressure can thus be decomposed into 
its different components as 


Peg — Pdig + Pc + Prr 


(18) 


where the components are obtained explicitly 

df{K,d,A) 


PPdig — 


dd 


d 




rd/2 


PPrr 


I c{z)dz 

-d/2 

= n(0) - nb- 


- Co 


(19) 

( 20 ) 
( 21 ) 


The contribution PfUg arises from self-interactions of ran¬ 
dom charges on the surfaces of the two slabs with their 
(dielectric and/or salt) image charges. This contribution 
stems from the first term of the interaction free energy 
(fT4l) and is present irrespective of multivalent ions, be¬ 
ing typically comparable in range and strength with the 
vdW pressure. It can, however, scale differently with the 
inter-surface distance and can be repulsive or attractive 
depending on the sign of dielectric discontinuity param¬ 
eter, A (e.g., it is repulsive for A > 0, applicable to the 
cases studied in this paper), resulting thus in a rather di¬ 
verse behavior for the inter-surface interaction between 
two net-neutral slabs in the absence of multivalent ions 
[ 75 I [ 7 ^ . The two terms Pc and Pmon correspond to the 
disjoining or osmotic pressure contributions of multiva¬ 
lent and monovalent ions, respectively. The former stems 
from the second term of the interaction free energy da, 
while, as noted above, the latter is included heuristically 
through Eq. (1^ . 

The total pressure acting on the slabs is finally given 

by 


P = Peg + PvdW- (22) 

This quantity is plotted in Fig. [3^ as a function of the 
rescaled inter-surface distance for fixed x = 2, cq = 1 mM 
and four different values of uq = 20, 25, 30 and 35 mM, 
which correspond to kIb = 0.35, 0.38, 0.42 and 0.45, re¬ 
spectively (colored dashed curves). For these same pa¬ 
rameter values, we also show the vdW pressure in Fig. 




(b) 


FIG. 3. (a) Rescaled total interaction pressure as a function 
of the rescaled distance between the randomly charged inner 
surfaces of two net-neutral dielectric slabs for fixed x = 2, 
Co = 1 mM and no = 20, 25, 30 and 35 mM as indicated on 
the graph (colored dashed curves). The corresponding vdW 
pressure for these parameter values closely overlap (black 
dashed curves). Panel (b) shows the total pressure (Ptot) 
along with its electrostatic {Pea) and vdW {Pvdw) compo¬ 
nents for Co = 1 mM and no = 20 mM. The light-colored 
dashed curves show Pea for no = 25, 30 and 35 mM (from 
right to left). 


[3K (black dashed curves), but these latter curves closely 
overlap and are not discernible at the implied resolution. 
This clearly indicates that the dependence of the total 
pressure on the salt screening parameter enters mainly 
through the electrostatic contribution (see below). In 
fact, as one can see from the figure, the total pressure, 
P, is less attractive than the standard vdW pressure at 
large separations, but, as uq decreases, the total pressure 
becomes more strongly attractive as compared with the 
vdW contribution. We show the electrostatic and vdW 
components of the total pressure separately in Fig. [ 8)3 
for Co = 1 mM and no = 20 mM {kIb = 0.35). The 
electrostatic contribution, P^g, shows a non-monotonic 
behavior: It is repulsive and decays to zero at large sepa¬ 
rations, while it becomes strongly attractive and diverges 
at small separations, showing thus a weakly repulsive 
maximum at intermediate inter-surface separations. The 
light-colored dashed curves show the same quantity for 

















(a) 




(b) 


(c) 


(e) 





0.01 


- 0.01 - 


d/U 

(d) 



FIG. 4. Rescaled electrostatic pressure as a function of the rescaled distance between the randomly charged inner surfaces of 
two net-neutral dielectric slabs for fixed kl^ = 0.35 (co = 1 mM and no = 20 mM. Panels (a) and (b) show the results for 
weakly and strongly disordered cases in a dielectrically homogeneous system (A = 0), while panels (c) and (d) show the same 
for A = 0.95. The results for the disorder-free case (black solid curve) and the vdW pressure (black dotted curves) are shown 
for comparison. Panels (e) and (f) show the contributing components to the electrostatic pressure, Pes, in the examples of a 
weakly disordered system with x=0.5 and a strongly disordered system with x=2, respectively. The other parameters are as 
in Panels (c) and (d). For comparison, we also show the contribution of mono- and multivalent ions in the disorder-free case, 
Pmon and Pc ■ Insets show the rescaled pressure components over a smaller range of values. 


larger values of no = 25,30 and 35 mM (from right to 
left). Clearly, the extent of the repulsion and the non¬ 
monotonic behavior of the electrostatic pressure become 
more pronounced as no is increased. Non-monotonic in¬ 
teraction pressures due to image-induced ion depletion 
have also been found in the absence of surface charge 
disorder but at relati vely large bulk salt concentrations 
(e.g., above 250mM) (T^ . 

Before proceeding with a more detailed analysis of the 
electrostatic contribution, we emphasize here that Pes 
and Pvdw show clearly different qualitative behaviors as 
one can see from Fig. [3 |d. The standard DLVO de¬ 
scription of colloidal interactions ii in terms of the 
vdW interaction of neutral surfaces is therefore insuf¬ 
ficient when the surfaces are neutral only on the aver¬ 
age and otherwise carry random positive and negative 
charges, and/or when the system contains also an asym¬ 
metric Coulomb fluid. Our results show that only a small 
degree of charge randomness with, e.g., a surface charge 
variance oi g = 0.04 nm“^ (corresponding to y = 2), 
will be enough to generate a sizable deviation from the 
standard vdW prediction. 


C. Interplay between charge disorder and 
multivalent ions 

In order to gain further insight into the intriguing role 
of disorder-induced effects in the presence of multivalent 
ions, we examine the behavior of the effective electro¬ 
static pressure, Pes, and its components in more detail. 

First, we consider two strictly neutral (disorder-free) 
surfaces, in which case the disorder effects vanish and 
the electrostatic inter-surface pressure, Pes, is given only 
by the osmotic pressures of multivalent and monovalent 
ions. The results are shown by black solid curves in 
Figs. [Ha, b) and Figs. He, d) for a dielectrically homoge¬ 
neous (A = 0) and a dielectrically inhomogeneous system 
(A = 0.95), respectively. For the latter case, the con¬ 
tributing components of the electrostatic pressure, i.e., 
P° and Pg^^, are plotted in Figs. He, f) for comparison. 
The vdW pressure is shown by black dotted curves in 
panels a to d. In all cases, we have taken cq = I mM and 
rig = 20 mM. 

It turns out that the electrostatic pressure in a 
disorder-free system (y = 0), which is given by Pes = 
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P^on + j is negative (attractive) for all inter-surface 
separations and tends to zero at large separations be¬ 
cause of salt screening effects. This means that the neg¬ 
ative bulk pressure is stronger in magnitude than the slit 
pressure for both cases of A = 0 and 0.95. This is be¬ 
cause of the image-induced depletion of ions generated 
by salt and/or dielectric-image charges that are of the 
same sign when the surfaces are immersed in a medium 
of larger dielectric constant as assumed here. In the limit 
d —>■ 0, Pes in an inhomogeneous system reduces to the 
bulk pressure, Pj, = {ni, + Co)k^T, due to complete de¬ 
pletion of ions from the slit under the strong repulsive 
force of dielectric-image charges, while in a dielectrically 
homogeneous system, complete depletion is not achieved; 
hence, Pes is less attractive in this latter case. On the 
other hand, for the given parameter values, the vdW 
pressure is more (less) attractive than the electrostatic 
pressure, Pes, for A = 0.95 (A = 0) as one can see from 
the figures. 

Next, we discuss how the introduction of charge disor¬ 
der on the inner surfaces of the net-neutral dielectrics af¬ 
fects the effective electrostatic interaction between them. 
To this end, we vary the disorder strength parame¬ 
ter from X = 0 up to X = 5 (corresponding to g = 
0 — 0.1 nm“^) and divide the results into two categories 
of “weakly” and “strongly” disordered cases. 

In the weak disorder regime (typically X ^ l)j f^e 
electrostatic interaction pressure, Pes, becomes gradul- 
lay more repulsive or positive (especially at small sepa¬ 
rations) as compared to that of strictly neutral surfaces 
when X is increased up x = 1 (see colored dashed curves 
in Figs. [4^ and c for A = 0 and A = 0.95, respec¬ 
tively). The repulsive interaction pressure is stronger 
for larger values of the dielectric discontinuity parame¬ 
ter, A (compare panels a and c). The most remarkable 
feature of our results is that Pes in this regime adopts 
a non-monotonic behavior: For very small Xj it devel¬ 
ops a point of zero electrostatic pressure followed by a 
shallow attractive minimum at intermediate separations 
(see, e.g., x = 0.25, red dashed curve, in panel a, or 
X = 0.25 and 0.5, red and blue dashed curves in panel c). 
The depth of this minimum decreases and the interac¬ 
tion profile instead develops a repulsive (positive) hump 
at large values of x (e.g., X = !)■ These behaviors are in 
clear contrast to what is observed in the case of strictly 
neutral (disorder-free) surfaces (black solid curves). 

In order to elucidate the origin of this difference, we 
compare different components of Pes in the case of an 
inhomogeneous system with A = 0.95 and x = 0.5 in 
Fig. 0)2. For comparison, the light-blue solid curve and 
the black solid curve (which is seen more clearly in the 
inset) show the osmotic contribution of monovalent and 
multivalent ions, P'^on S'lid P°, in the case of strictly 
neutral surfaces (x = 0). Both contributions are attrac¬ 
tive (negative), while the former is the dominant one as 
expected since we have assumed that multivalent ions en¬ 
ter in small bulk concentrations (here, cq = I mM). By 
introducing a weak charge randomness (x = 0.5 for the 


relevant curves in Fig. [4j;), the negative pressure of both 
multivalent {Pe, red dashed curves, inset) and monova¬ 
lent ions {Pmon, blue dashed curve, main set) decrease in 
magnitude and, thus, become less attractive as compared 
with the disorder-free case. It is important to note that 
the distribution of both monovalent and multivalent ions 
and, thus, their respective osmotic pressure, Pmon and 
Pc, are affected by the surface charge disorder through 
the disorder-induced, single-ion potential, which is the 
second term in Eqs. or m- The decrease in the 
magnitude of attractive pressure components Pmon and 
Pc in the weak disorder regime is because more mono- 
and multivalent ions are attracted into the slit from the 
bulk solution in the presence of surface charge disorder 
(see Fig.[2]and Fig. Obelow). Nevertheless, at very small 
separations, all ions are again totally depleted due to the 
repulsive forces of image charges and Pmon and Pc re¬ 
duce to the same bulk values as P^^^ and do when 
d ^ 0, i.e., —rihk^T and —c^k^T, respectively. The 
crucial role, however, is played by the repulsive disorder 
self-interaction component, Pdis (green dashed curve), 
whose effect is amplified by increasing the charge dis¬ 
order strength and/or the interfacial dielectric mismatch 
as follows from Eq. ([19]). It is the interplay between this 
repulsive disorder-induced pressure (which dominates at 
small to intermediate separations) and the attractive os¬ 
motic pressure of ions (which dominate at larger sepa¬ 
rations) that gives rise to the non-monotonic behavior 
mentioned above. For sufficiently small x, pressure 
of multivalent ions is typically too small (as compared 
to the other two components) to change this qualitative 
behavior. 

This picture changes drastically in the strong disor¬ 
der regime (typically X ^ l)i where the pressure due to 
multivalent ions becomes a key factor. In Fig. [Hd and 
d, we plot Pes for larger disorder strength parameters, 
X = 1.5 up to 5 for A = 0 and A = 0.95, respectively. 
In this regime, we find a reverse trend caused by the 
surface charge disorder: The repulsive hump at interme¬ 
diate separations now diminishes and eventually disap¬ 
pears when X is increased to larger values; one thus finds 
a highly attractive (negative) interaction pressure with 
a range that can be much larger than that of the vdW 
interaction pressure (black dotted curves), or the mere 
image-induced, ion-depletion pressure (black solid curves 
with X = 0), at sufficiently large disorder strengths. 

The different components contributing to the electro¬ 
static pressure for this case are shown in Fig.|4f for x = 2 
with other parameters being the same as Figs. [H: and 
d. As seen, while the disorder self-interaction contribu¬ 
tion, Pdis, has become only slightly more repulsive (green 
dashed curve), the other two components. Pc and Pmon, 
show significant changes as compared to their counter¬ 
parts in the weak disorder regime in Fig.jU^. Specifically, 
the attractive contribution due to the osmotic pressure 
of multivalent ions (red dashed curve) becomes much 
stronger than what we find in the weakly disordered or 
disorder-free cases, while the contribution from mono- 
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d/lg 

FIG. 5. Average total number of multivalent ions in the slit, 
N, per rescaled surface area, S = S/l%, is plotted versus the 
rescaled inter-surface separation, d/lB, for a system similar to 
those considered in Figs. ID: and d. 

valent ions now becomes repulsive (positive) in contrast 
to its attractive (negative) behavior in these latter cases 
(compare panels e and f). The monovalent contribution 
becomes positive because of an increase in the number 
of monovalent ions in the slit that creates a larger os¬ 
motic pressure. The average number of multivalent ions 
also increases in the slit as x increases; however, these 
ions correlate strongly with the disorder charges on the 
bounding surfaces of the slabs (through the second term 
in Eq. (fl^ which contributes directly to Pc, Eq. (l20l) l 
and, as a result, give rise to an effective attraction be¬ 
tween the net-neutral slabs through the pressure compo¬ 
nent Pc- This effect is in clear contrast with the image- 
induced, ion-depletion mechanism, in which a decrease 
in the number of mono- and multivalent ions in the slit 
gives rise to attractive osmotic pressures on the slabs 
similarly from both types of ions. This behavior origi¬ 
nates from the (singular) attractive, single-ion potential, 
which is created by surface charge disorder. It is thus one 
of the most fundamental aspects of the coupling between 
surface charge disorder and mobile multivalent ions that 
follows from our results. 

The remarks in the preceding discussions on the at¬ 
traction (depletion) of multivalent ions to (from) the slit 
region can be corroborated by considering the average to¬ 
tal number of ions in the slit. For multivalent ions, this 
quantity can be calculated from N = S It 

shows a very different behavior as a function of the inter¬ 
surface separation in the weak disorder and strong dis¬ 
order regimes as one can see in Fig. [5] (we have used the 
same parameters as in Figs.lH: and d). Although, at fixed 
inter-surface distance, always more multivalent ions are 
attracted to the slit by increasing the disorder strength, 
the response of multivalent ions to the decrease in slab 
separation is quite different for different disorder regimes: 
For disorder-free and weakly disordered surfaces, mul¬ 
tivalent counterions are quickly depleted from the slit 
by decreasing the slab separations, while, for sufficiently 
large disorder strengths, these ions (as well as monova- 




FIG. 6. (a) Rescaled total inter-surface pressure as a function 
of the rescaled distance between the randomly charged inner 
surfaces of two net-neutral dielectric slabs for Hxed A = 0.95, 
X = 2, no = 100 mM and co = 1, 5 and 10 mM as indicated 
on the graph (from top to bottom). The corresponding elec¬ 
trostatic and vdW components of the total pressure, i.e., Pes 
and P-udw, are shown in rescaled units as well (Pes is shown 
by the light-colored dashed curves and Pvdw by the black one, 
which closely overlap). Panel (b) shows Ptot along with P^dw, 
Pea and the three components contributing to the latter, i.e., 
Pdia, Pc and Pmon for fixed Co = 10 mM and other parameter 
values as in (a). 


lent ions that are not shown here) are more strongly at¬ 
tracted to the slit from the bulk solution due to stronger 
attraction experienced from the randomly charged, inner 
surfaces of the slabs. In the intermediate regime of disor¬ 
der strengths, one find a non-monotonic behavior upon 
deceasing the inter-surface separation, with the average 
total number of multivalent ions first decreasing and then 
increasing at very small separations. One should, how¬ 
ever, note that in dielectrically inhomogeneous systems 
such as the one considered in Fig. [5] with A = 0.95, 
dielectric-image repulsions eventually win and, thus, even 
the curves for highly disordered systems eventually turn 
downward and we find iV —>■ 0 for d —>■ 0 or, equivalently, 
Pes —t —{nb + co)kBT for all cases in Figs.[4l: and d (these 
limiting behaviors occur at very small separations that 
are not physically meaningful, e.g., below d ~ 1 A, and, 
for the sake of presentation, they have not been shown 
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in the plots). 

The crossover from weak to strong disorder regimes 
can be effected also by decreasing the salt screening pa¬ 
rameter through a decrease in the bulk concentration of 
monovalent ions, hq (and/or increasing the bulk concen¬ 
tration of multivalent ions, cq). For instance, at suffi¬ 
ciently large uq such as no = 100 mM, the electrostatic 
interaction pressure, Pes, at the disorder strength value 
of X = 2 and multivalent salt concentration of cq = 1 mM 
shows a weak-disorder behavior (see Fig. [^), in contrast 
to what we found for y = 2 at lower salt concentration of 
rig = 20 mM in Fig. |4ji. Even though increasing cg up to 
cg = 10 mM enhances the attractive pressure as shown in 
the figure, the behavior of the effective interaction still 
remains in the weak-disorder regime as can be verified 
from the pressure components in Fig. [6 |d (e.g., compare 
Pmon and Pc with those in Fig. [djs). This is clearly be¬ 
cause of strong salt screening effects at large ng, in which 
case the strong-disorder behavior can be achieved only by 
taking even larger values of y. 


V. CONCLUSION AND DISCUSSION 

In this work, we have studied the effective interactions 
between net-neutral dielectric slabs that carry quenched 
(fixed) random charge distributions on their apposing 
surfaces while they are immersed in an asymmetric aque¬ 
ous electrolyte (Coulomb fluid) containing mobile mono¬ 
valent and multivalent ions. The effective interaction be¬ 
tween quenched, random charge distributions have been 
considered in a series of recent works [s^, Issl - I^ . IH- 
1^. ItsI - I^ . 1^ . and, while the role of Debye screening 
due to a weakly coupled (monovalent) salt solution on 
these interactions have been analyzed, the role of multi¬ 
valent ions that couple stron gly with surface charges has 
received much less attention (621. WA . 

The interaction between (electro-) neutral dielectric 
slabs is traditionally described in terms of vdW or 
Casimir-type forces [l|, that are always attractive be¬ 
tween identical bodies in vacuum or in a polarizable 
medium such as water or an aqueous salt solution. In 
the DLVO context , the vdW attraction is counter¬ 
acted by the mean-field electrostatic repulsion between 
like-charged surfaces, providing thus a mechanism for the 
stability of colloidal dispersions. In the case of strictly 
neutral surfaces, however, the electrostatic inter-surface 
repulsion is obviously absent but, in addition to the vdW 
interaction, there are image-induced, ion-depletion forces 
that are attractive as well and arise because the mobile 
solution ions are depleted from the vicinity of the di¬ 
electric interfaces as a consequence of the repulsion from 
their image charges (these image charges are of the same 
sign when the surfaces are immersed in a medium of 
larger dielectric constant as compared with that of the 
slabs). This effect is intrinsically non-mean-field and oc¬ 
curs because of the discrete nature of ions, neglected in 
the collective mean-field description based on the stan¬ 


dard Poisson-Boltzmann theory Image-induced, 

ion-depletion effects for strictly neutral dielectric surfaces 
have been studied extensively in recent years (see, e.g.. 
Refs, [sl-l^ and references therein) and the role of mul¬ 
tivalent ions, in particular, has been considered in Refs. 

mm- 

In the present work, we add a new feature to the cur¬ 
rent understanding of the effective interaction between 
neutral surfaces and consider the situation, in which the 
surfaces are neutral only on the average but otherwise 
carry a quenched random distribution of positive and 
negative charges. Indeed, heterogeneously or randomly 
charged surfaces are commonplace in soft matter and 
have attracted a lot of attention in recent years (see, 
e.g.. Refs, [dsl - l^ l. We have thus shown that the pres¬ 
ence of quenched surface charge disorder in conjunction 
with mobile multivalent ions in the ionic solution leads 
to remarkable and significant qualitative changes in the 
standard picture commonly accepted for the interaction 
between neutral bodies. 


The charge disorder has several manifestations in the 
present context: First, the self-interactions between the 
quenched random charges on the inner surfaces of the 
dielectric slabs and their image charges lead to a repul¬ 
sive, short-ranged interaction pressure, Pdis, that tends 
to counteract the attractive vdW interaction pressure, 
Pvdw as discussed thoroughly elsewhere [751 - 1^ . The 
surface charge disorder, on the other hand, modifies the 
single-ion interaction potentials in a way that creates a 
singular attractive potential acting between ions and the 
dielectric boundaries 1^, 0, 


84| . This effect works 


against the ion-image repulsions that tend to deplete ions 
from the slit region between the slabs. Thus, as a re¬ 
sult of this attractive disorder-induced potential, more 
mono- and multivalent ions are accumulated in the slit 
when the confining dielectric boundaries are randomly 
charged. However, we find qualitative differences in the 
ways this excess attraction affects the osmotic pressure 
of ions. The osmotic pressure due to monovalent ions, 
P-mom can become repulsive which, in fact, agrees with 
the standard mean-field picture since monovalent ions are 
only weakly coupled to the bounding surfaces and create 
a larger entropic pressure upon further accumulation in 
the slit. On the other hand, the osmotic pressure due to 
multivalent ions. Pc, can become even more attractivel 
This is in clear contrast with the image-induced, ion- 
depletion mechanism, in which a decrease in the number 
of mono- and multivalent ions in the slit gives rise to 
(stronger) attractive osmotic pressures on the slabs, sim¬ 
ilarly for both types of ions. 


This rich behavior stems from combined effects of sur¬ 
face charge disorder and the presence of mobile multi¬ 
valent ions, with the latter creating strong inter-surface 
attractions upon further accumulation in the slit. As 
a result, the effective total interaction pressure between 
randomly charged, net-neutral dielectric surfaces can dif¬ 
fer qualitatively from what one expects based on the 
standard vdW and/or image-induced, ion-depletion in- 
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teractions for strictly neutral surfaces, where they predict 
purely attractive interactions (unless at relatively large 
bulk salt concentrations, e.g., above 250mM, as noted in 
Ref. [13 )• In particular, the net effect from the compet¬ 
ing electrostatic components due to disorder, mono- and 
multivalent ions results in a distinct, disorder-induced 
non-monotonic behavior. 


In the weak disorder regime, the net electrostatic pres¬ 
sure becomes gradually more repulsive (positive) at in¬ 
termediate to large separations as the disorder coupling 
parameter (or equivalently the disorder variance, g) on 
the bounding surfaces of the slabs is increased from zero 
(strictly neutral or disorder-free system), exhibiting first 
a shallow minimum and then a pronounced repulsive 
hump at intermediate separations. As the disorder cou¬ 
pling parameter is further increased into the strong dis¬ 
order regime (typically beyond x = ~ 1), the 

behavior of the net electrostatic pressure is reversed; the 
repulsive hump is diminished and one finds a strongly 
attractive (negative) interaction pressure with a range 
and magnitude larger than that of the vdW or the 
image-induced, ion-depletion interaction pressures be¬ 
tween strictly neutral surfaces. The crossover from weak 
to strong disorder regimes can be effected also by decreas¬ 
ing the salt screening parameter (through decreasing the 
bulk concentration of monovalent ions) and/or increasing 
the bulk concentration of multivalent ions. The two dif¬ 
ferent regimes of the disorder effects, the strong and the 
weak disorder regimes, bear also some resemblance to the 
strong-weak coupling dichotomy that exists for the net- 
charged su rface s, dependent on the electrostatic coupling 
parameter 103l| . It is no small feat to be able to partition 
the behavior of this complicated system by two dimen¬ 
sionless coupling parameters that effectively govern its 
salient characteristics. 


Finally, we note the qualitative differences between the 
ion-mediated interactions we find between net-neutral 
surfaces in this work and those we reported between ran¬ 
domly charged surfaces carrying a finite mean surface 
charge density, cr 0, in Refs. [H, [13] • One of the 
key differences is that the repulsion between mean sur¬ 
face charges on the inner surfaces of the slabs, Po-, in 
the latter case contributes a dominant repulsive pressure 
to the net electrostatic pressure, which is stronger, both 
in magnitude and range, than Pdis and Pvdw- On the 
other hand, when the surfaces carry a finite mean sur¬ 
face charge density, the osmotic pressure of monovalent 


ions, Pmom turns out to be always repulsive and com¬ 
parable in magnitude and range to P^', these repulsive 
pressure components completely mask the short-range in¬ 
teraction pressures, Pdis and Pvdw, whose roles in com¬ 
petition with the attractive pressure due to multivalent 
ions, Pc, are brought up only in the case of net-neutral 
surfaces. Notably, we find that Pmon hr the case of sur¬ 
faces carrying a finite mean surface charge density hardly 
responds to changes in the disorder strength, while, in 
the case of net-neutral surfaces, it can change from a 
purely attractive component (due to image-induced, ion- 
depletion effects) for disorder-free and weakly disordered 
surfaces to a repulsive one (due to disorder-driven ion ac¬ 
cumulation in the slit) for strongly disordered surfaces. 
Another difference between the two cases is that multi¬ 
valent (counter-) ions in the case of surfaces carrying a 
finite mean surface charge density correlate strongly also 
with the mean charge densities of the confining bound¬ 
aries of the slit, in line with the strong-coupling paradigm 
for electrostatics of charged surfaces known to generate 
stro ng l ike-charged surface attractions (see, e.g., Refs. 
[9^104 |1. Such strong-coupling effects are clearly absent 
in the present case with net-neutral surfaces. 

In summary, our analysis provides new insight into the 
intricate role of surface charge disorder in the context of 
net-neutral surfaces and its fingerprints on the effective 
ion-mediated interactions between them. These interac¬ 
tions are predicted to occur with a range and magnitude 
comparable to or, for strongly disordered systems, much 
larger than the vdW and/or image-induced, ion-depletion 
effects considered previously only in the case of strictly 
neutral surfaces. Our predictions will thus be amenable 
to numerical and/or experimental verihcation. 
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